Forward-Backward rapidity correlations at all rapidities 
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The study of correlations and fluctuations can give an important contribution to the understanding of the behaviour 
of dense matter in high energy heavy ion collisions. See, for recent experimental developments [1,2,3,4] and for recent 
Qh' theoretical studies [5,6,7,8]. 

We consider forward-backward (FB) rapidity correlations in the framework of a two-step scenario [9,6,7,8] with 
formation of longitudinal sources, as it happens in the glasma [9] or in the dual string model [10]. Even in the case 
OC of particle emission from the sources being fully uncorrelated, long distance E-B correlations are unavoidable due to 
CI fluctuations in the colour or the number of sources. In the situation of multiple elementary collisions, as it occurs in 
nucleus-nucleus collisions, and hadron-hadron collision, as well, such fluctuations are naturally expected [11,12]. 
The basic formalism for the development of the two-step scenario is as follows. In a high energy nucleus-nucleus, 
Cln' or hadron-nucleus, or hadron-hadron collision the interaction occurs via v elementary interactions with formation of 
JL. , sources: longitudinal fields or strings. In the dual string model- which we shall use as reference, but nothing essential 
• depends on that - the number N s of formed strings is directly related to the number of elementary collisions: 

N s = 2v . (1) 

Next, the strings emit particles as a local, short range correlated in rapidity, phenomenon. If P(n Fl n B ) is the 
■ probability of emitting n F particles in the backward one, we have 

G\ ■ 

P(n F ,n B ) = ^2ip{v) ^2 J\p(n>iF,ni B )S(n F -^2n iF )S(nB -^2n iB ) (2) 

V {n iF ,n iB } i i i 

o ; 

where <p{y) is the probability of occurence of v elementary collisions and p(ni F ,ni B ) the local probability of nip 
particles in the forward window and ms particles in the backward one from the i-th string. 
^ , We make now the usual assumption that particle emission from the string is not long range correlated such that, if 
a, large enough, rapidity gap exists between the forward and the backward window, then 



p(n iF ,n lB ) = p(n iF ) p(n iB ) (3) 

and, for the average of ni F rii B , 

n lF n iB = nTJ" . (4) 

However, local short range correlations in the forward (and backward) window exist - due, for instance, to reson- 
nances - and characterized by the local variance 

d2 iF = n \ F - n*F 2 ■ ( 5 ) 
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If the strings are non-interacting and are all identical and symmetrical with respect to central rapidity - which 
requires symmetrical AA collisions and flat rapidity distribution - then in any rapidity window one has all the N s 
strings contributing which means that in the summations J2 i and in the products Y\ { , i takes the values, see (I), 



i = l,2,...,2i/. (6) 

However, in general, (6) is not correct, as the number of strings contributing in a given rapidity window is, in 
general smaller than 2v. There are two reasons for that. First, strings may fuse and percolate such that their number 
decreases [13,14]. Second, some strings are smaller than others - that is the reason why the particle density in rapidity, 
dn/dy, is not flat - which implies an additional rapidity dependent reduction. These reduction factors are naturally 
functions of the window rapidities, yp and y B , respectively. 

Instead of the limit 2v in (6) we thus write for the number of contributing strings 



F : 2v — ► 2vr{y F ) , (7) 

B : 2v — > 2vr{y B ) , (8) 

where r(jjp) and r(y B ) are the reduction factors, with < r(yp), r(jjB) < 1. As one imediately sees, we do not have 
to make estimates on r(yp) and r(y B ) has they appear via the measurable multiplicities (rip) and (ns). 
With (3), (5), (7) and (8) in (2) we obtain: 

(n F ) = 2<v> r(y F )n F , (9) 

(n B ) = 2<v> r{y B )n B , (10) 

r>2 / \ / w \ ( n F)(n B ) nn 

D FB = {n F n B ) - (n F ) (n B ) = — , (11) 



and 



Dip = (n%) (n F f = M! + f- (up) , (12) 



K np 

where 1 /K is the normalized fluctuation of the v distribution: 



(1/2) - („) 

Usually, one defines the F — B correlation parameter b as 



, n ; ; , , 2 • (is) 



^4, (14) 



FF 



such that, from (9), (10), (11), (12) and (14), we obtain 



ft . (15) 

1+ K dp V ' 

~ (n F ) n F 

With symmetrical situation, AA collisions, and y B = — yp, (n B ) = (np), we obtain the well known result [15,7]: 



1+ * ^£ 

(hf) np 



(16) 



It is clear from (15) that there exists the possibility of b being larger than 1. This is particularly true if one neglects 
short ramge fluctuations: dp = 0, and 
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FIG. 1: a) STAR central data, 0-0.6%, Au - Au data at RHIC, and our curves (full lines) for D 2 FB (11) and D 2 FB (12). 
b)Predictions for LHC, Pb — Pb central (dashed lines). 



b ~ (n B ) I (n F ) . 

In general it is assumed that local emission follows a Poisson distribution, such that d F — n F and 



(17) 



b = 



(n B ) I (n F ) 
1 + K/ (n F ) 



(18) 



One should keep in mind that (16) and (17) are quite general, in the sense that they do not require symmetric AA 
collisions or symmetric windows. 

The framework of the two-step scenario with long range correlations due to extended sources in rapidity and local 
short range correlations, is also present in the colour Glass Condensate (CGC) approach [6,8]. The result obtained 
has the expected structure of (15) or (18): b = A[l + B\~ l . 

Let us test the validity of (18) for b and (11) and (12) for D 2 FB and D FF , respectively, in the symmetrical situation, 
AA collisions and y B = —y F ((n F ) = {n B ))- The STAR collaboration at RHIC has presented new results for D 2 FB 
and D FF in central collisions at mid rapidity, making use of windows of width Srj = 0.2 and a distance between their 
centres of Ar/(0.2 < A77 < 1.8) [3]. In (11) and (12) K has fixed at the value of 88 1 . 

In Fig. la) we compare (11) and (12) with mid rapidity data for central Au — Au at y's = 200GeV and show our 
expectation in forward rapidity region. In Fig. lb) we show our prediction for Pb — Pb at LHC, y/s — 5500GcV. In 
Fig. 2 we compare b, (18) with RHIC data and we show as well our predictions for LHC. Note that our prediction 
qualitatively agrees with the CGC model [6]. In our plots we have used for highest centrality particle density the 
parameterization of [16], 



1 dn 
A~chj 



„-(l- Q) y 



+ 1 



(19) 



where Y is the beam rapidity. The parameters in (19) where ajusted to fit the PHOBOS/RHIC dn/drj data at 
V^s = 200GeV [17]: A = 0.26, a = 0.31, 5 = 0.75. 



1 As the STAR value for D pB went up by a factor of around 4 [3], the value of K used in [7], goes down by the same amount. 
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FIG. 2: STAR data and our curve (full line) for b, (18). Predictions for LHC, Pb — Pb central (dashed lines). 



Regarding the parameter K, (13), which is defined at partonic/string level, we would like to mention that, in the 
case of A A collisions, K is related to a similar parameter K n , defined at nucleon level and which can be estimated 
in a Glauber calculation [18]. As K n depends mostly on the nucleus-nucleon cross-section and a nn slowly varies with 
energy, one expects K n to be fairly independent of energy. On the other hand, repeating the arguments of [12], which 
essentially says that the normalized fluctuation at the parton level should be larger than the same quantity at nucleon 
level, we obtain: 



K<K n ~ const. (20) 

The first relation in (20) is satisfied [18]. 

Having in mind (20) and the fact that K is a KNO (Koba, Nielsen, Olesen) invariant quantity [19], we have assumed. 



K = const., (21) 

independent of energy, for our predictions at LHC, in Figs.l and 2. Note that K, similarly to K n , is expected to 
increase with the number of participating nucleons. 

In the general case of asymmetrical AB collisions and asymmetrical windows one can use the parameterizations of 
[20]. 

There is a simple way of testing the model for forward-backward correlations in the general case (asymmetric 
distributions and asymmetric windows). Let us fix ys, the centre of the backward window, and move uf, the centre 
of the forward window, with xjf > Vb- We can rewrite (18) in the form 



1 + K'x 



(22) 



with 



< x < oo, and 



x = (n B ) I (n F ) , 



(23) 
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FIG. 3: The correlation parameter 6 as a function of x = (ub) / («f) ■ The function (21) is universal, only depending on the 
value of K' = Kj (hb) ■ In the figure K' = 1. 



K' = K/{n B ) . (24) 

Independentely of the shape of the rapidity density distribution, of the centrality, of the energy, of the rapidity size of 
the windows, the functional form of the dependence of b on x is given by (22). Of course K' is a number that varies 
with centrality (number of participating nucleons), energy and window size. However, the shape of the b dependence 
on x is always given by (22). In Fig. 3 we show the dependence of b on x (in the case of K' = 1). 

Finally, instead of studying forward-backward correlations, backward-forward correlations can as well be studied. 
One has simply to make the changes B — > F and F — > B in our formulae. 
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